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Abstract 

There are processes in nature that resemble a true force but arise due to the 
minimization of the local energy. The most well-known case is the exchange in- 
teraction that leads to magnetic order in some materials. We discovered a new 
similar process occurring in connection with an electronic phase separation tran- 
sition that leads to charge inhomogeneity in cuprate superconductors. The min- 
imization of the local free energy, described here by the Cahn-Hilliard diffusion 
equation, drives the charges into regions of low and high densities. This mo- 
tion leads to an effective potential with two-fold effect: creation of tiny isolated 
regions or micrograins, and two-body attraction, which promotes local or intra- 
grain superconducting pairing. Consequently, as in granular superconductors, the 
superconducting transition appears in two steps. First, with local intra-grain su- 
perconducting amplitudes and, at lower temperature, the superconducting phase 
or resistivity transition is attained by intergrain Josephson coupling. We show 
here that this approach reproduces the main features of the cuprates phase dia- 
gram, gives a clear interpretation to the pseudogap phase and yields the position 
dependent local density of states gap A(r) measured by tunnelling experiments. 
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1. Introduction 

Almost 25 years of intense research on copper-oxide-based high-temperature 
superconductors has revealed many interesting and non-conventional results. It is 
impossible to explained all these data in detail by a single theory because some 
features may be sample dependent, but there are some general properties believed 
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to be common to all cuprate superconductors like the pseudogap phase and the 
d-wave superconducting amplitude symmetry[[li0,[3l]. On materials which allow 
surface studies and thin films, mostly Bi, La and Y-based cuprates, there are well 
measured properties: the larger gap at the leading edge of the Fermi surface or 
antinodal (along the Cu — O bonds) direction ((±;r,0) and (0, ±7r))|3], with the 
consequent nodal Fermi arcs above the critical temperature Tcdsila], and the spa- 
tially dependent local density of states gap A(r) measured by atomically resolved 
spectroscopy such as scanning tunneling microscopy (STM)lf7l. [sl. lol. [lO . H, 121. 

We show here that the fundamental superconducting interaction and some of 
these properties can be interpreted evoking an electronic phase separation (EPS) 
and its consequent temperature evolution of the charge inhomogeneity. There 
are many evidences that the charge distribution in the Cu02 planes of the high 
temperature superconductors (HTSC) is microscopically inho mog eneous f, 1 3l 1 141 . 
Several different experiments like neutron diffraction 11 1 SL 1 1 6l. 1 1 711 . muon spin re- 
laxation (/i5i?) ||l8l |l91jnuclear quadrupole resonance (NQR) and nuclear mag- 
netic resonance NMR ||20L l2ll. l22|] have detected varying local electronic densities. 
The spatial variations of the electronic gap amplitude /^scij) at a nanometer length 
scaleflBB [10LlUJ,ll2|] may also be connected with the charge inhomogeneities. 
The origin of this electronic disorder is still not clear; it may be from the quenched 
disorder introduced by the dopant atoms[21, 23], or it may be due to competing 
orders in the CuOi planes ||24i] . Probably, the most convincing evidence of an EPS 
transition came from the increase of the local doping difference with decreasing 
temperature in the NQR experiment of Sing er et al[20]. An EPS has been also 
used to interpret transport properties 11251 12611 and to provide a clear interpretation 
to the non-vanishing magnetic susceptibility measured above the resistivity tran- 
sition rc[|27i|281]. More recently a phase separation in La2Cu04^y was analyzed 
after different times of the annealing processes ll23ll 

The cuprates phase diagram has a crossover temperature or upper pseudo- 



gap detected b 
doped region y. 



j_everal transport experiments that merges with Tc in the over- 
i.e., for large dopant average level p. It is our fundamental 
assumption, and the starting point of our theory, that such crossover line is re- 
lated to the EPS transition temperature Tps{p). Phase separation is a very gen- 
eral phenomenon in which a structurally and chemically homogeneous system 
shows instability toward a disordered composition[29]. However, the problem is 
to describe and to follow quantitatively the time and temperature evolution of the 
charge separation process. We have shown that an appropriate way to do this is 
through the general Cahn-Hilliard (CH) theory HEEQ- 

This approach uses an order parameter that is the difference between the tem- 
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perature dependent local doping concentration p{i,T) and the average doping 
level p, i.e., u{i,T) = {p{i,T) — p)/p. Then the local Ginzburg-Landau (GL) 
free energy functional is, 

f{iJ)=^-e'\Wu[iJ)\' + VGL{u{iJ)). (1) 

Where the potential Vgl{uJ,T) = -A^{T)u^ /A+ ...,A^{T) = a{Tps{p)- 
T), a and 5 = 1 are temperature independentparameters. e gives the size of the 
grain boundaries among two distinct phases f 3 ll. 13211 . 



The CH equation can be derived[29] from a continuity equation of the local 
free energy /(z, T), dtu = —V.J, with the current J = M'V{df/du), where M is 
the mobility or the charge transport coefficient, normally incorporated in the time 
step intervals. Therefore, 

1^ = -MV^e^V^u + VcdpJJ)). (2) 
We have already made a detailed study of the CH differential equation by finite 



difference methods p l|l which yields the density profile u{p,i, T) in a 105 x 105 



array as function of the time steps, up to the stabilization of the local densities [l32l 



33Ll34j,l35|,|3^. Here we study the EPS profile as the parameter A changes from 
A = near T = TpsioA= \ close ioT = OK. 

In Fig.([T]) we show a typical density profile after a long time evolution (A,5 = 
1 and average doping p = 0.16) with a high segregation level as it is seen by 
the local densities histogram (in the inset). The colour map shows regions with 
p{i) ~ (red) and p{i) ^ 2p (blue). The lack of direct measurements does not let 
us know whether all cuprates have such high level of charge inhomogeneity and it 
is possible that it varies according to a specific family of compounds. 



2. The EPS and the Fundamental Interaction 

During the EPS process the holes move preferable along the Cu — O bonds 
or nearest neighbor hopping according to the measured dispersion relations [l37|l . 
Here we show by the local free energy calculations that this hole motion can be 
regarded as if it was produced by an effective hole-hole interaction as has been 
proposed before, like, for instance, by the work of Trugman[38]. To see this we 
perform numerical simulations of the local potential energy VGL{u{i, T)) with the 
solutions u{i,T) of the CH equation. The results are shown in Fig.© and they 
reveal two important and different effects: /- it divides the system in many tiny 
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Figure 1: (color online) The order parameter u{i) = {p{i) — p)/p simulation of the inhomoge- 
neous local charge density on a 100 x 100 sites. In the inset, the local densities p{i) histogram of 
the optimal compound (p = 0.16) indicating the tendency to segregate into grains with p{i) w 
and p{i) k, 0.32. As the temperature T goes down below the phase separation temperature Tps the 
segregation increases. 



potential wells as it is shown on the free energy map of Fig.©. The inset shows 
also the free energy along the white straight line on 42 sites (~ 16QA) showing 
these potential wells and the barriers between the low and high density regions. 
We define the height of these inter-grain or grain-boundary potential V^b since it 
gives a granular structure to the system. In this way each small grain of low or high 
density is a small bounded region by Vg^. Each of these potential wells form local 
single-particle bound states. These bound states are seen experimentall y by the 
local density of states (LDOS) derived from STM experiments flBBEillllES] • 



Vg}j{T) is function of the temperature and increases when the temperature goes 
down below Tps as it is demonstrated by the two insets in Fig.© taken from 
simulations at two different temperatures. 

//- Second, in the process of minimizing the local free energy, the holes move to 
cluster themselves in a similar fashion as if they attract themselves, forming the 
hole rich and hole poor regions. 

This hole-hole effective attraction is schematically illustrated in Fig.® where 
a) represents a homogeneous system with p = 0.25 (one hole at each four sites) 
and b) the motion toward clusters formation of low (p, = 0) and high densities 
{pi = 2p = 0.50). This hole movement can be regarded as originated from an ef- 
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Figure 2: (color online) On the top panel, the potential V{p,i,T) planar map simulation. The 
values on 42 sites at the white line are shown in the top inset to demonstrate the potential wells 
with average barrier ss Vgt, where single-particle bound states appear. The lower panel is a 3- 
dimensional section of the above map, where the potential barriers are more clearly visible, at a 
higher temperature and smaller Vgt,, as shown in the inset. 



fective two-body attraction. Conceptually, this is similar to the spin-spin exchange 
interaction, that arises from the Pauli principle and the minimization of the local 
electronic energy, that produces ferromagnetic order. 

We made simulations using VQi{p,i,T) at low temperatures up to Tps{p), 
where the disorder disappears. The low temperatures values are parameterized 
to yield the low temperature gaps at different average doping p from the STM 
measurements of the series of Bi2212 LDOSf?]. In this way we obtain a numer- 
ical estimate of Vgb{p,T) that reproduces the STM results of this entire family. 
Namely 

Vgt{p.T)=V{p)V{T) = (-0.7 + 2.4p)(l-(r/rp5)i-^), (3) 
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Figure 3: Schematic representation of the cluster and pair formation, a) represents a homogeneous 
system with p = 0.25 (one hole at each four sites) and the CH phase separation produces hopping 
preferably along the Cu — bond toward clusters formation of low (/?,■ = 0) and high densities 
ipi = 0.50)(b). At low temperatures, d-wave superconducting pairs may be formed (c). 



where the values are in eV , V{p) is linear and vanishes aX p 0.32 following 
the approximate behavior of Tps{p) from the 7b (p) curve plotted in the phase 
diagram of Timusk and Statt[l]. The power to the 1.5 temperature dependence is 
taken from the work of Cahn and Hilliard[30]. 



3. BdG-CH Combined Calculations 

We now perform self consistent-calculations with the Bogoliubov-deGennes 
(BdG) theory in the system. There are some possibilities to obtain a solution 
with a disordered density and local dependent superconducting gaps, for instance: 
Ghosal et al[[39|l introduced a varying impurity potential in the chemical potential 
to assure the local changes in the hole density. Other approach is to introduce 
modulations to the superconducting pair interaction due to the out of plane dopant 
atoms [40]. Here we introduce a third method, different but it contains the spirit of 
the work of Ghosal: we keep the local disordered density solution derived from 
the CH simulation fixed at all points and determine self-consistently the respective 
local chemical potential and d-wave superconducting amplitudes. In this way 
we capture the phase separation solution like that shown in Fig.©. To calculate 
the superconducting amplitudes and single-particle excitations and bound states 
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we used the effective hole attraction Vgt, in the form of a local two-body nearest 
neighbor attraction. Although the potential is constant in space, the different local 
densities yield different superconducting amplitudes. Starting with an extended 
Hubbard Hamiltonian with nearest neighbor hopping t = 0. 15eV and next nearest 
neighbor hopping of ti jt = -0.64, close to the ARPES value of ti jt = -O- SOpTII 



the BdG mean-field equations are written in terms of the BdG matrix OZ |33 



34 



35L136|] 



These equations, defined in detail in Refs.l 32l 34], are solved self-consistently 



for En > together with the eigenvectors {umV„) and d-wave pairing amplitudes 
Ad(x,-) = -^£[m„(x,)v*(x/ + 5)+v*(x/)m„(x,- 

^ n 

+5)].anh^, (4) 

and the input inhomogeneous hole density is given by 

p(x,) = 1-2£[|m„(x,0|V. + |v„(x,0|'(1-/„)], (5) 

n 

and converges self-consistently toaN xN square (here we made calculations with 
N = 2S — 42) in the CH density map of Fig.Q. /„ is the Fermi function. Vgb is 
maintained fixed at each temperature T and for a given compound with average 
doping p. 

The BdG-CH combined calculations yield larger superconducting amplitudes 
(/,/), T) in the high density grains and smaller amplitudes at regions with low 
densities. The variations on Arf(z,p) are usually ±12% around the average value 
A2^'(r,p) = Aj{T,i,p)/N. However, due to the mean-field approach, all the 
amplitudes vanishes at the same temperature T*{p). We believe that a more rigor- 
ous treatment would yield that larger gaps vanish at larger temperatures, but there 
are many interesting consequences even within this simple approach. Due to the 
functional form of Vgb{p) in Eq.© the A'^^{T,p) decreases systematically with p 
as it is shown in Fig.©. As the temperature decreases below T*{p) the grains 
become superconductors but the grain boundary potential barrier Vgt, prevents the 
current to flow freely. Thus, the electronic grain structure of cuprate supercon- 
ductors at low temperatures can be regard as formed by numerous 5,- — / — Sj 
junctions. 

Consequently, we apply the theory of granular superconductors ['47| to these 
tiny electronic grains, and as an approximation, we use the Josephson coupling 
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expression to an S — / — S junction given by [|41ll . 



A^^'(r,p) is the average of the calculated BdG superconducting gaps A^i{i,T) on 
the entire x A'^(A^ = 28 — 42) square. The Rn is the normal resistance of a given 
compound, which we take as proportional to the pab measurements [42] on the 
complete series of Laj-pSvpCuOi- The values of A'^^{T,p) as function of p are 
plotted in the top panel of Fig.©. In the low panel we plot the two sides of 
Eq© and the intersections yield one of our main result: the Tc{p) dome shape in 
agreement with the Bi2212 series. 

The zero resistivity transition takes place when the Josephson coupling Ej 
among these tiny grains is sufficiently large to overcome thermal fluctuations [4311 . 
i.e., Ej{p, T ~ Tc) ^ ksTcip), what leads to phase locking and long range phase 
coherence. Consequently the superconducting transition in cuprates occurs in two 
steps, similar to a superconducting material embedded in a non superconducting 



matrix [|43ll : First, as the temperature goes down, by the appearing of intragrain su- 
perconductivity (pseudogap phase) and by Josephson coupling with phase locking 
{superconducting phase). We emphasize that the reason for the dome shape form 
of Tc{p), with the optimum value near p = 0.16, is due to the fact that /S!^J{T,p) 
decreases and {l/Rn) increases with p. 

4. The STM Results and Interpretation 

One of the most important experimental result without a widely accepted ex- 
planation is the spatial dependenj: energy gap s A(r) measured by Scanning Tun 



neling Microscopy (STM)|S i, Isi OJ, 0, [120. We show here that this behavior 
can be reproduced by the microscopic granular theory developed above. We per- 
form calculations on the non-uniform charge system but the usual local density 
of states (LDOS)[44] at different places x,- are not so sharp and difficult to get 
their real value. To improve the determination of the peaks positions, although 
the measured LDOS of cuprates are not symmetric, we deal with the following 
symmetric LDOS, 

Ni{E)=Y, [ |m„(x,)|2 + |v„(x,)P]x 

n 

[f'„{E-E„)+fl{E+E„)]. (7) 
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Figure 4: In the top panel we plot the average Af/(r,/?) and in the inset the schematic single 
particle energies that appear in the LDOS. In the low panel, the Josephson coupling among super- 
conducting grains Ej{p,T) for some selected doping values as function of T and ksT. The curves 
intersections give the dome shape Tc{p), as plotted in the inset. 



/„ is the Fermi function, the prime is the derivative with respect to the argument, 
and M„, v„ and are respectively the eigenvectors and positiv e eig envalues (quasi- 
particles exciting energy) of the BdG matrix equation OZ l33l. l34i |35L 13611 . 

To study the effects of the disordered density, we examine the ratio LDOS{Vgb ^ 
0)/LDOS{Vgb = 0). LDOS{Vgh = 0) contains the inhomogeneous charge distri- 
bution and LDOS{Vgh 7^ 0), the unnormalized local density of states, vanishes 
around the Fermi energy due to the superconducting gap and the single particle 
levels in the grains. This LDOS ratio yields well-defined peaks and converges to 
the unity at large bias, in close agreement to similar approach of the LDOS ratios 
calculated from STM measurements [9]. 

These features are illustrated in Fig.® for an overdoped p = 0.20 compound 
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near a grain boundary at a representative average doping hole point. AtT = 40K, 
Vgh = 0.240eV and for large applied potential difference (V > 0. leV) both LDOS, 
LDOS{Vgh 7^ 0) and LDOS{Vgb = 0) converge to the same values, but around zero 
bias {V ~ 0) there is a large spectral weight suppression in LDOS{Vgh only. This is 
due to the single-particle bound states and to the superconducting pair formations 
and both effects are destroyed as the temperature is raised, according to the lower 
panel of Fig.©. To distinguish clearly the effect of the superconducting gap in 
the DOS of Fig.® we draw arrows at the values. Notice that it appears only 
as an anomaly in the LDOS curve with two well-defined peaks, as it was reported 
by some STM data[|lOl. 1 1111 . Consequently, the gap calculated directly from the 
LDOS peaks is identified with the with the pseudogap Apq{p). This also implies 
that STM may not measure the local superconducting gap but, in some cases, the 
local gap from the single particle bound states formed in the isolated islands. 




-0.06 -0.04 -0.02 0.02 0.04 0.06 
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Figure 5: The real LDOS with V,,t = 0.240eV, the zero pair potential LDOS Vgt = and the 
normalized LDOS on an overdoped p — 0.20 sample for y^,t = 0.240eV. The zero pair potential 
LDOS Vgb = is also shown. The arrows show the values of the local superconducting gap 
A,/ ^ l.lmeV, smaller than the LDOS gap (w 24meV). 



In Fig.® we plot the calculated low temperature LDOS at four representative 
locations with different densities of the p = 0. 16 optimum compound. We see that 
the peaks vary from 35 — 65meV in accordance with the data of McElroy et al[0]. 
This overall agreement comes from the choice of the parameter in the potential 
Vgb (Eql3]), however the local variations in nanometer scale comes from the charge 
inhomogeneity calculated from the CH approach and also reproduces the it local 
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variations of the STM LD0S[7]. The temperature evolution in the low panel of 
Fig.® shows that the single-particle or LDOS gap closes near T* = 1 15K, more 
than 20K above Tc. 




Figure 6: In the top panel, the LDOS at four different regions with local doping n,. At the low 
panel the temperature evolution of the LDOS just to the closing of the gap at T « 1 14K, above 
Tc w 92K. The arrows show the average A"i^'{T — 4QK) ~ l5meV. In low panel we show the LDOS 
temperature evolution at one region (pi = 0.25). The kink at low bias in the low temperature curve 
has also been observed by many STM data and in our calculations it is due to the 

superconducting gap Aj(r = 40K) = 2lmeV. 



The calculations in the density disordered system yield that the low tempera- 
ture superconducting gaps Arf(z, p = 0. 16) vary from 7 — 23meV, depending on the 
local density, that is, they are much smaller than the gaps derived from the LDOS 
peaks shown in Fig.®. Consequently we conclude that the measured LDOS gaps 
are due to the single-particle bound states gaps ApQ{i,p, T) that leads to the fol- 
lowing scenario to the cuprate superconductors: The pseudogap phase derived 
from our calculations is formed due to the granular behavior by the single-particle 
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bound states gaps ApQ{i,p,T) in the isolated local doping dependent regions or 
grains. Like in granular superconductors, there are also the intragrain supercon- 
ducting amplitude A^(z,p, T). At lower temperatures the superconducting phase 
also contains both gaps but has phase coherent through the Josephson coupling 
among the grains. At the overdoped region T*{p) approaches Tc{p) while their 
difference increases in the underdoped region. These two calculated curves as 
function of p are shown in the phase diagram of Fig.©. 

We have also plotted the EPS line Tps{p) used in our calculations where stats 
the phase separation process and which we assumed to be near the observed 
anomaly calledupper-pseudogap Timusk and Statt[[ll]. As mentioned in the in- 
troduction, this transition may be due to the presence of the AF order which has 
lower free energy ('2?] (intrinsic origin) or due to out of plane arrangements of 
dopant atoms Dili |22|,|23fl (extrinsic origin). 




0.05 0.1 0.15 0.2 0.25 0.3 
P 

Figure 7: The calculated phase diagram of cuprate superconductors as derived from the EPS 
transition Tps{p) and the formation of tiny grains. The onset of single particle bound states and 
superconducting amplitude formation occurs at T*{p). The dome shape Tc{p) curve is due to 
Josephson coupling among the grains. These three lines are in agreement with the experimental 
results. 



5. Conclusions 

The main ideas of our work come from the application of the CH theory to 
model the highly non-uniform charge distribution in cuprates. This approach de- 
scribes the formation of local isolated regions of low and high densities formed 
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at the many free energy potential minima. These non-uniform regions act as at- 
traction centers and reduces the charges kinetic energy what favours the Cooper 
pair formation. We show that an effective hole-hole attraction Vgh{p) appears 
from the Ginzburg-Landau free energy and the Cahn-Hilliard simulations. These 
simulations reveal also the microscopic granular structure that leads to the charge 
confinement and the single-particle bound states (as seen in Fig. (2)). With the 
values of Vgh{p) matching the low temperature Bi2212 average LDOS gap values 
of McElroy et alR and the assumption of long range order by Josephson cou- 
pling among the different regions we reproduce the entire phase diagram shown 
in Fig.©. The dome shape form of Tc{p) is a consequence of the different behav- 
ior of 1 /Rn{p) which increases, and A^^'(p) which decreases with p. 

In summary: the LDOS main peaks are likely to be due to the intragrain 
single-particle bound states Apo and the local superconducting gap A^ produces 
in general only a small anomaly at low temperature and low bias as observed 
by many STM datali, [HI [M US]- A^ and A^g have completely different na- 
ture but they have the same origin, namely, the inhomogeneous potential Vgi,. 
They are both present in the pseudogap and superconducting phases but they have 
distinct strengths. Their differences were verified in the presence of strong ap- 
plied magnetic fields and also by their distinct temperature behaviour in tunneling 
experiments [45]. The potential barriers Vgb prevents phase coherence and, as in 
a granular superconductor, the resistivity transition occurs due to Josephson cou- 
pling among the intragrain superconducting regions (Eq©. 
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